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Abstract 

In these lecture notes I describe some of the main theoretical ideas 
emerged to explain the aging dynamics. This is meant to be a very short 
introduction to aging dynamics and no previous knowledge is assumed. I 
will go through simple examples that allow one to grasp the main results 
and predictions. 
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1 Contents 



As the title makes it clear, the aim of these lectures is not to present an exhaustive 
and complete description of the state of the art in aging dynamics (even less an 
exhaustive reference list). But rather, its purpose is to present in a nutshell 
the main theoretical ideas emerged to explain the aging dynamics. For a more 
complete and deep presentation there are already very good reviews and lectures 
in the literaturePU Oj. These references are also useful to have an introduction 
to the experimental results on aging (in particular the ones on memory and 
rejuvenation effects) which will not be discussed in the following. 

The contents of these lectures are: 

• Introduction to aging and off-equilibrium thermal relaxation 

• A simple example: domain growth in the Ising model 

• Mean-Field Theory of Aging 

• Activated dynamic and scaling 

• The Trap model 

• Outstanding issues 
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2 Introduction 



Off-equilibrium dynamics is a very wide and interesting subject. Many works 
have been devoted to it recently new interesting theoretical concepts as well as 
fascinating experimental results have been discovered. Yet, it is certainly fair 
to say that the overall comprehension is still not complete and certainly many 
new discoveries (as well as much work) is ahead. There are roughly two different 
types of off-equilibrium dynamics: driven off-equilibrium dynamics and thermal 
off-equilibrium relaxation. In the first case the system is driven out of equilib- 
rium by an external stationary force (it may be a shear for a liquid, a voltage for 
an electron system, etc. . . ) and kept in a stationary off-equilibrium state. The 
subject of these lectures is the latter type of off-equilibrium dynamics: thermal 
off-equilibrium relaxation. In this case the system, at equilibrium at time to, 
is brought out of equilibrium by changing a control parameter C (for example 
the magnetic field for a magnetic system, or the temperature for a glass-forming 
liquid, etc. . . ). Subsequently the system, which is in general in contact with a 
thermal bath, starts to evolve towards the new equilibrium state corresponding 
to the new value of C . However, this may take a long time which actually can 
become infinite in the thermodynamic limit or just for practical purposes. 
From a theoretical point of view this is a clear example in which "more is differ- 
ent". Indeed a system with few degrees of freedom coupled to a thermal bath is 
going always to equilibrate in a finite time 1 (which in the worst case scenario and 
for the majority of systems scales with the number N of the degrees of freedom 
as e NK where K is a positive constant). However, when the number degrees of 
freedom is very large this equilibration time may become infinite for all practical 
purpose. As in standard equilibrium statistical mechanics the way to deal with 
this situation is study the system taking the thermodynamic limit from the be- 
ginning thus focusing on the evolution of the system on timescales which may be 
extremely large but do not diverge with the system size. 

Let us give just an example before introducing the aging dynamics. Consider a 
three dimensional Ising ferromagnet at a temperature T larger than the critical 
temperature T c . If one quenches the system from T to T' > T c the ferromagnet 
will relax (i.e. its energy, its local correlation functions, etc. . . , will relax) on a 
finite timescale to the new equilibrium state at temperature T'. However, it has 
been proven that for T' < T c the largest relaxation time scales as expKN 2 ^ 3 , 
hence, in the thermodynamic limit, the system remains out of equilibrium forever 
(in the next section we are going to understand in detail the properties of the off- 
equilibrium dynamics in this case). Note that this never-ending off-equilibrium 
regime is quite general whenever the quench crosses a phase transition. Let us 
consider another example: the measurement of thermo-remanent magnetization 

1 We neglect the cases in which the configuration spaces of the system with few degrees of 
freedom is broken in disconnected pieces. 
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Figure 1: Thermo-remanent magnetization M, normalized by the field-cooled 
value Mf c , vs. t(s) (log 10 scale) for a spin-glass (Ag : Mn2s% sample), at T = 
9K = 0.87T C . The sample has been cooled in a 0.1 Oe field from above T g = 
10. 4K to 9K; after waiting t w , the field has been cut at t — 0, and the decaying 
magnetization recorded (from [2]). 

in a spin glass. In this case the system is cooled in a small field from above 
T c to a temperature below, it then waits in the field for a time t w after which 
the field is switched off. The subsequent evolution of the magnetization is very 
peculiar and an example is reported in Fig. ^ After a very rapid fall-off of the 
magnetization there is a very slow decreasing. A very important feature of these 
evolutions is that the longer is t w the slower is the relaxation. This behavior is 
called aging because the system evolution depends on its age (t w ). This means 
that the typical relaxation timescale for the system is not fixed a priori but is 
evolving and is fixed by the age of the system itself: older systems relax more 
slowly. 

In order to understand better this point let us again focus on a system with a 
finite number of degrees of freedom [which therefore has a finite relaxation time] . 
Generically one expects a similar aging behavior for not too large t w . However 
when the typical relaxation time scale at time t w becomes of the order of the equi- 
libration time all the aging effects are interrupted and the system is characterized 
by an equilibrium dynamics. This in particularly means that one time quanti- 
ties like the magnetization, the energy, etc. equal their equilibrium values and, 
hence, they are time independent; whereas for example the two-time correlation 
and response functions are invariant under translation of time, i.e. they depend 
just on the time difference. Finally typical equilibrium relations like fluctuation 
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dissipation relations between correlation and response are verified. 
Thus, it is clear that aging and the never-ending relaxation discussed before is a 
collective phenomenon with (at least one) growing correlation length 2 . Indeed 
one of the theoretical approaches to understand and explain aging dynamics is 
based on determining what is the growing correlation length, how it grows with 
time and assuming scaling with respect to this length. Another approach is based 
on the analytical solution of infinite dimensional (or also infinite number of com- 
ponents) models. In this case an interpretation in terms of a growing lengthscale 
is at a first sight indirect. Instead the standard physical interpretation that has 
been developed during the last ten years is in terms of evolution within the energy 
landscape. 

In the following we shall analyze the aging dynamics in the Ising model. This 
is a good example to start with because it is rather simple and allows one to 
introduce both approaches which are correct and complementary (in this case at 
least!). In the following sections I will present more in detail both approaches 
focusing in particular on their application to disordered systems. Finally I will 
also briefly present the trap model, that although certainly not a microscopic 
model, is used often in phenomenological studies and appears as a coarse grained 
description for some microscopic models. 

3 A simple example: domain growth in the Ising 
model 

In the following I shall focus on the aging dynamics of the three dimensional 
Ising model due to a quench from the high disordered temperature phase to a 
temperature at which, in equilibrium, the system is ferromagnetically ordered. 
This example is particularly enlightening because it can be analyzed in a sim- 
ple approximated way (of course much more sophisticated analyzes have been 
performed). Furthermore all the different theoretical ways to tackle the aging 
dynamics can be introduced and their physical content can be easily understood. 
The first thing to do in order to understand what is going on is to look at the 
results of a numerical simulations. In Fig. |21 there are two snapshots of a 2d 
cut of a 3d ferromagnetic Ising model after a quench from the high to the low 

2 Why? Well, it is the usual sloppy argument. If all the correlation lengths are bounded 
during the relaxation then one can roughly divide the system in independent sub-systems 
for which the relaxation time is bounded by a certain timescale ts- Thus the aging has to be 
interrupted and the system has to relax on a timescales not larger than ts which contradicts the 
hypothesis of never reaching equilibrium on any finite timescales. In real systems this argument 
can be used to see that one should have a growing (although not diverging) correlation length 
even for systems with interrupted aging. Finally, a word of caution, this argument might 
have problems for system with long range interactions, mean- field models and quench to zero 
temperature. 
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Figure 2: Two snapshots of a 2d cut of a 3d ferromagnetic Ising model evolving 
with a Glauber dynamics after a quench at time t — from high temperature 
to T < T c . Black color is used for minus spins and white for plus spins. On 
the left the configuration reached after t w = 10 3 MonteCarlo steps. On the right 
the configuration reached after t w = 10 5 MonteCarlo steps. Statistically the two 
configurations looks the same after a length rescaling. 



temperature phase (black color is used for minus spins and white for plus spins). 
This figure clearly shows that the systems is trying to separate minus and plus 
regions in order to gain bulk energy and decrease interface energy: the system is 
separated in domains that grow in time. 

In order to analyze this off-equilibrium dynamic regime it is convenient to set 
up a continuum description in terms of the coarse-grained magnetization density 
4>{x, t) which obeys a Langevin equation: 



d t (p(x, t) = 



dF 



(1) 



where r/(x, t) is a white noise with a variance determined by the temperature T 

< n(pc, t)r)(x!, t') >= 2T5(x - x!)8(t - t') 
and F is the standard Ginzburg-Landau free energy 



(2) 



The other important information that one has to provide is the initial condition 
that we are going to take of the form: 



< 0(x, O)0(x', 0) >= A5(x - x') 



(3) 



Since we are focusing on a system that is in equilibrium at a very high temperature 
at t — 0, the correlation length is extremely small and, at our coarse-grained 
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level, the magnetization density correlations are just given by Q. At t = + the 
bath temperature is instantaneously switched to T < T c , where T c is the critical 
temperature of the Ising model. Hence, the system starts to evolve trying to 
thermalize without never succeeding on any finite time. 

In this case the temperature T does not play a very important role, as long as 
it is smaller than T c . Therefore we shall put it equal to zero and we will discuss 
later what happens for non zero temperature. 
As a consequence we have to solve the differential equations: 

«9 t 0(x, t) = (A + m 2 )0(x, t) - g<f>{x, t) 3 (4) 

with the initial conditions (JBJ). These eqs. correspond to a pure gradient descent 
in the energy landscape defined by F. There are two terms in F: the bulk-energy 
term — ^(j) 2 + |0 4 that is minimized by (p = ±^Jm 2 /g and the (V0) 2 which is 
minimized by flat configurations. As a consequence the system during its gradient 
descent will try to decrease F and this will lead to fatter and flatter regions with 



a = ±ym 2 /g as seen in Fig. |21 

Equations (j3J) are somewhat simpler than the original stochastic ones but they 
are still too complicated to be analyzed exactly. Therefore in the following we 
are going to make use of an approximation a la Hartree that can be justified as 
a 1/N expansion. It consists simply in replacing the nonlinear term <?0(x, t) 3 by 
g < 0(x, t) 2 > 0(x, t) where < • > means the average over the initial condition 3 . 
Now the equations read simply: 

d t <f>{x,t) = [A + a{t)]<f>{x,t) (5) 

where a(t) = m 2 — g < 0(x, t) 2 > has to be determined self-consistent ly (note 
that since the initial average is translation invariant < 0(x, t) 2 > does not depend 
on x). Eq. (JSJ) can be immediately integrated in Fourier space: 



(k, t) = 0(k, 0) exp (^-kh + f* a(t')dt' 



and the self-consistent equation on a(t) reads: 

a{t) =m 2 - gA J exp (-2k 2 t + 2 J* a{t')dt^j 

where eq. (J3J) has been used to eliminate the initial condition 4 . Since we are 
interested only to large times and we know that the systems is trying to decrease 

3 If one generalizes the field </>(x, t) to an n-components field tfi a (x,t) and F to 
J fix |(V0) 2 — I] j-$ 2 + 4n($ 2 ) 2 then the equation 10} becomes exact in the infinite N limit. 

The reason is that in eq. the term g<j> 3 becomes g4> a (f> 2 /N and in the infinite N limit the 
term <j) 2 /N does not fluctuate and is equal to its average. 

4 Some readers can find that this expression is a bit singular at t = 0. In this case one has 
to remember that there is an underlying lattice and there is a cut-off in the integration over k 
corresponding to the inverse of the lattice spacing. 
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its local energy — + |0 4 as well the elastic ones ^(V0) 2 we expect (and it 
can indeed be shown) that at large times 

< 2 (x, t) >= A J exp (-2k 2 t + 2 jf* a(t')dt') 

converges to a finite value (actually equal to m 2 /g). But this is possible only if 

2 /* a(t'W ^ 3/2 log(*/* ) *o = {Am 2 / g) 2/3 /8tc 

JO 

at large times, i.e. a(t) ~ ^ for t » 1. Now that we know the expression of a 
we can compute the space time correlation function at large times: 

This expression is very interesting and deserves different remarks. First, for any 
finite and fixed (x — x'), it becomes at large times simply a function of t'/t. Thus, 
we have found that the correlation between two points, as far as they could be 
in space and time, converges to the value m 2 / g when £',£—>• oo provided that 
their distance in space and time is kept fixed when the large time limit is taken, 
i.e. if one looks to a finite region on a finite time this region will be typically 
on one of the two state ±^m 2 /g. On the other hand, for fixed but very large 
time t, t' the magnetization density correlation always fall to zero over distance 
|x — x'| oc y/t + f. This means that, at time t, the typical size of the regions in the 
states ±\Jm 2 /g is L(t) oc y/t. Furthermore, for fixed |x — x'|, the magnetization 
density correlation always falls to zero over a time separation t — t' which is larger 
than t' (we assume t' < t). This means that the time it takes to the system to 
decorrelate from its configuration at time t' is of the order of t' itself, i.e. the 
age of the system is the characteristic timescale for the dynamical evolution: the 
older is the system, the slower is its dynamics. 

Thus we recover what discussed in the introduction and plotted in Fig. El the 
systems is decomposed in domains of size L(t) (that is the growing correlation I 
talked about in the introduction), for which we have found a growth law L(t) oc 
y/t. In the interior of each domain the system is in one of the two state ±^Jm 2 /g. 
However, if one waits a time of the order of the age of the system a point x has 
been swept by different domain walls and thus the magnetic correlation has been 
lost. 

In this approximate treatment of this simple model we find the essential ingre- 
dient of the coarsening/scaling picture of the aging dynamics. On long time-scales 
the system is broken up in domains which grow and the dynamics is self-similar 
provided the length is appropriately rescaled by L(t). This length L{t) is natu- 
rally interpreted as the typical size of the domains and its rate of growth can be 
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found from general arguments In this case L(t) oc \fi but the growth rate 
may be different for different systems j^j. For example, if we consider the Ising 
model but with a dynamics that conserves the magnetization, for example the 
Kawasaki dynamics [3], then L(t) oc t 1 ^ 3 . However, because of the self-similarity 
of the dynamics, one generically expects that the correlation function can be 
written in a scaling form 



This is indeed the case in our approximate treatment and we refer to E] for a 
more general discussion. 

Another way to think about aging and glassy dynamics is the one coming from 
the solution of mean-field disordered systems very much based on the energy 
landscape (or more precisely free energy but here it does not matter because 
T = 0). Although we postpone a more technical discussion to the next section, 
in the following we shall introduce the main ideas PQ 0] . Let us start with some 
general remarks. 

The equation (J3J) consists in a gradient descent in the energy landscape from a 
high temperature configuration toward the ground state configurations (that are 
never reached on finite time scales). During this gradient descent the velocity 
of the system decreases and tends to zero in the infinite time limit. This rather 
natural results can be checked noticing that the square of the absolute value of 
the velocity in the energy landscape v = — VF = — ^|^y equals Since 
on general grounds one expects that one time quantities, in particular F, tends 
toward a well defined value in the infinite time limit (indeed one can check that 
this is the case within our Hartree approximation) we obtains that and the 
velocity has to tend to zero in the infinite time limit. 

Another important quantity other than the velocity that is useful to track the 
system evolution is the the energy Hessian evaluated at the configuration reached 
at time t, in particular its density of eigenvalues. Since the system is trying to 



reach the completely stable states ±ym 2 /g we expect that the fraction of negative 
eigenvalues becomes smaller and smaller at large times and vanishes at t = oo. 
Let us check this assumption within our approximate treatment in which the 
Hessian is diagonal in Fourier space: 






Computing the density of states is an easy job and one finds: 
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Figure 3: Pt(A) as a function of A for t = 1,2,5,10 (from left to right). For 
simplicity iV is put equal to one. 



where N is a normalization factor 5 . Since m 2 — g < 0(x, t) 2 >= a(t) ~ 3/4t we 
find that the density of states of the Hessian has at its left edge a square root 
form starting at a negative value \ m in — — 3/4t which indeed approaches slowly 
zero at infinite times, see Fig El Finally, the last piece of information that we 
need is that the velocity always points toward directions corresponding to small 
eigenvalues of the Hessian. This can be checked noticing that when the energy 
decreases slower than exponentially to its infinite time value (a quite general 
behavior in aging system that is true in our case): 



d?F 
dt 2 
dF 
dt 



VF-H • VF 

- 2 ^vf? ,0 



(8) 



when t — > oo. Hence, decomposing this identity along the eigenvectors of H and 
calling the components of VF on the eigenvector fi we find: 







which shows that indeed the velocity always points toward directions correspond- 
ing to small (vanishing) eigenvalues of the Hessian. 

Putting all these pieces of information together, we obtain an energy landscape 
picture of aging in which the system makes a gradient descent in the energy 



5 Note that since there is an underlying lattice there is also a cutoff in the integration over 
k. Therefore N is cut-off and lattice dependent. Also p(X) is in reality cut-off and lattice 
dependent except at its left edge where is given by eq. Q. 
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landscape and at long times it moves slowly and slowly (the absolute value of the 
velocity is going to zero) along channels with many uphill directions and a low 
density of downhill directions. The system velocity points towards the directions 
over which the energy landscape is flatter. The longer is the time the smaller 
is the fraction of downhill directions and they are more and more flat. Thus it 
takes longer to the system to move and go downhill. This never-ending descent 
in the energy landscape dominated by an increasing flatness of the landscape is 
the main picture behind aging obtained from the exact solution of the dynamics 
of different mean- field models pQ. 

Finally, let me comment on how things change if the final temperature after the 
quench is not zero (but still less than T c ). In this case inside the domains the 
system has a pseudo-equilibrium dynamics that takes place on a finite timescale. 
For example, the local spin-spin correlation function will be approximatively a 
function of t — t' for t — t' « t' (and t' » 1) that converges to a plateau mf on 
a finite timescale where is the local magnetization. The aging regime shows 
up only on timescales such that t — t' oc t' . It corresponds to the exit from the 
plateau and can still be described in the same way discussed above. 

4 Disordered systems 

In the following I will discuss in some detail the theoretical approaches introduced 
to explain the aging dynamics. I will focus mainly on disordered systems for 
historical reasons and simplicity and I will comment in the conclusion on the 
generalization to non-disordered glassy systems. 

4.1 Mean-Field Theory of Aging 

The study of the equilibrium and off-equilibrium dynamics of mean-field disor- 
dered systems revealed that a very rich and interesting behavior can be found 
within the mean-field approximation. Generically, it has been found that these 
systems displays a dynamical transition at a temperature at which the re- 
laxation time diverges. A quench to a temperature less than leads to an off- 
equilibrium behavior and aging that persist forever if the thermodynamic limit 
has been taken from the beginning. The way in which the dynamical transi- 
tion affects the equilibrium dynamics for T > is discussed in the Lectures 
of A. Cavagna |5] and D. R. Reichman jHj and I will not discuss it. I would 
like just to point out that two quite different behaviors may take place. In 
mean-field disordered systems having a spin-glass transition, like for example the 
Sherrington-Kirkpatrick model, the dynamical transition temperature coincides 
with the temperature at which the equilibrium glass transition takes place. In- 
stead, for systems like the p-spin disordered model, which have an equilibrium 
transition conjectured to be related to the structural glass transition, the dy- 
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namical transition temperature is larger than the equilibrium one. Furthermore, 
contrary to the first case, the transition has a (partially) first order character 
since a plateau strictly different from zero emerges in the correlation functions at 
Td, see the left plot in FiglU 
Historically, the mean-field theory of aging has been developed analyzing the 




Figure 4: Left: Sketch, for a p-spin model, of the decay of the stationary corre- 
lations in the high T phase close to Td, T\ > T% > . . . (from left to right). Right: 
Sketch, for a p-spin model, of the decay of the aging correlations in the low T 
phase, at fixed T < Td, t[ < t' 2 < . . . (from left to right). 

off-equilibrium dynamics of completely connected models. An example is the 
p-spins model whose Hamiltonian is 



Kii <...<i n 



and the couplings Ji lt ...,i p are independent Gaussian variables with zero mean 
and variance p\/2N p ~ 1 . The dynamics of this systems can be fully analyzed 
if the degrees of freedom are variables Sj subjected a global spherical constraint 
Y<iLi s 1 — N. The usual strategy is to take a Langevin dynamics for these (spher- 
ical) models and write down the Schwinger-Dyson equations for the correlation 
C and response functions R: 



1 ^ < Si(t) > 



c(t,t') = -Y,<s l (t)s t (t')> R(t,f) = -Y: hj[fl) • 

where hi(t') is a local magnetic time-dependent field coupled to and the overline 
means the average over disorder. A great simplification due to the completely 
connected character of these models is that the self-energy can be written as a 
simple polynomial function of correlation and response. More specifically, the 
general expression one finds is pQ: 

r+oo r+oo 

/ dt 1 RQ 1 (t,t 1 )C(t 1 ,t')-2TR{t',t) = / dt 1 Z c {t,t 1 )R(t',t 1 )+ (9) 
Jo Jo 
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r+oo 

Jo 

r+oo 

/ dhR^it^Rit^f) = 5(t-t')+ (10) 
Jo 

r+oo 

+ / dtJl R {tM)R{ti,t') 

Jo 

where Rq 1 is the inverse of the bare response function. For example in the 
case of the p-spin spherical model R^ (t,ti) = [d t + fi(t)]5(t — ti) where /i(t) 
is a spherical Lagrange multiplier determined enforcing the spherical constraint 
C(t,t) = 1. The self-energies are simple functions of the propagators and read: 

E c (t, t') = ]T c n [C(t, t')] n E R (t, = J2 nc n [C(t, tX^R^ f ) (11) 

n n 

and the coefficient c n depends on the particular model one is focusing on. For 
example in the case of the p-spin spherical model all the c n are zero except 
c p _i = p/2, i.e. c n = 

At this point it is easy to understand that it is possible to derive these type of 
equations not only within the exact solution of completely connected models but 
also within self-consistent approximate treatment of finite dimensional systems. 
Indeed in the case of Langevin dynamics one can set up a field theory in which 
the self-energy is written as a sum of diagrams constructed with the vertices of 
the theory and using the true propagator as line. The equation (jllj) is in general 
equivalent to retain only the first self-consistent diagram and can be often justi- 
fied in terms as a large N expansion in 1/N. This also makes it clear why this 
the type of eqs. introduced above can be found also for non-disordered systems. 
More complicated functional forms have been also partially studied in the liter- 
ature but I will not discuss it further. 

There are two methods for analyzing the eqs. (|9I1U|) : a complete solution can be 
obtained by integrating them numerically but much information can be gained 
already analytically. In the following I will just make a summary of the main 
findings. 

For T > correlations and responses converge at large times toward a time 
translation invariant regime (TTI), i.e. they depend only on the time difference 
t — t' and furthermore they are related by the fluctuation-dissipation relation 
R{t) = —C'(t)/T. Hence, the system is at equilibrium. Approaching Td the 
correlation function develops a plateau for discontinuous spin models and takes 
more and more time to attain the equilibrium TTI regime, see Fig. HJ Instead 
for T < Td the correlation and response never reach the equilibrium and a TTI 
regime. In simple cases, for example in the case of the p-spin spherical model, 
there are two time sectors: a pseudo-equilibrium one that corresponds to taking 
the limit t' — > oo with t — t' fixed. In this regime the correlation and response 
become indeed TTI and are related by the fluctuation dissipation relation (as for 
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the domain growth discussed in the previous section). The aging regime corre- 
sponds to t — > oo with t'/t fixed in which the correlation and response function 
become functions of t' / t. The correlation function equals the Edwards- Anderson 
parameter q^A for t'/t = 1 and approaches zero in the limit t'/t = (I take 
t' < t), see Fig |U Note that in more complicated case the t'/t is replaced by 
h(t')/h(t) with h a model dependent function. It may also happen, in particular 
for continuous models that there is an infinite sequence of time sectors, each one 
of them determined by a ratio hi(t') / hi(t) . 

Finally, a striking property of this off-equilibrium regime is that the correlation 
and response are still related by a fluctuation-dissipation like relation (FDR) in 
which, however, the bath temperature is replaced by an effective temperature 
T e ff. In the case of different time sectors there is an effective temperature for 
each one of them. The interpretation of T e ff as an effective temperature has been 
justified in different papers. Furthermore, it has also been shown that T e ff can 
be related in some cases to the Parisi function x(q) introduced in the study of 
the thermodynamics. This topic has generated an enormous amount of interest: 
many numerical simulations and some experiments have been performed in order 
to check the existence of this off-equilibrium relation between correlation and 
response. A very popular way to plot data to check for a generalization of FDR 
consists in plotting the integrated response — It w R(t,t')dt' as a function 

of the correlation C(t, t w ). If the fluctuation dissipation relation R = —C'/T e ff is 
verified (in the limit t w — > oo) one should find a straight line with slope —1/T e ff. 
In Fig. El we plot on the left a typical FDR plot for a system converging to equi- 
librium. In this case the curves, for different t w evolves toward a straight line. 
Instead for mean-field glassy systems off-equilibrium a rather different behavior 
has been found, see the right part of Fig. [5] The dot-dashed line corresponds 
to a case with an infinite number of timescales and T e ff. The continuous line 
with a breaking point corresponds to a systems, like the p-spin spherical model, 
with only one T e ff. The dotted line corresponds to the domain growth in the 
Ising model T e ff = oo. The other curves correspond to waiting times far from 
the asymptotic regime. A detailed discussion about T e ff can be found in [Tj. I 
will just make some other comments in the conclusion. 

Finally, let me discuss the physical interpretation of this aging dynamics. As 
anticipated in the previous section the interpretation that comes from many dif- 
ferent analytical and numerical studies is based on a never-ending descent in the 
(free-)energy landscape 6 . Indeed the aging after quench at zero temperature is 
very similar to his finite temperature counterpart. Instead a very different behav- 
ior would arise if aging was due to activated hopping over (free-) energy barriers. 
The interpretation is exactly the same outlined in the previous section on coars- 
ening except that now the eigenvectors of the Hessian change during the time 

6 In reality only at zero temperature this interpretation is fully justified because otherwise 
the dynamics cannot be mapped exactly into an evolution within the free-energy landscape. 
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Figure 5: The asymptotic behavior of the integrated linear response against 
the correlation in a parametric plot, for fixed waiting time and using t as a 
parameter. Left: behavior for a system relaxing to equilibrium. Right: behavior 
in a slowly relaxing system out of equilibrium. The dot-dashed line corresponds 
to a case with an infinite number of timescales h(t) and T e ff. The continuous line 
with a breaking point corresponds to a systems, like the p-spin spherical model, 
with only one T e ff. The dotted line corresponds to the domain growth in the 
Ising model. Finally the other curves correspond to waiting times far from the 
asymptotic regime. 



evolution. Another important new element with respect to domain growth in the 
Ising model is that for these mean-field glassy systems there are an exponential 
number (in the systems size) of thermodynamic states. It is generically believed 
and indeed it has been shown in simple cases that the system ends up, at long 
times, aging very close to the states that are marginally stable, i.e. that have a 
vanishing fraction of flat directions (at least after a very slow quench). Finally, in 
the case of discontinuous spin glasses, the effective temperature has been related 
to the complexity S(f), i.e. the log degeneracy of states with a given free-energy 
/, as 1/T e jf = dS/df where the derivative is evaluated at the free energy density 
of the marginally stable states. 

I will not continue further the discussion of the mean field theory of aging. I refer 
to more detailed review for deeper analysis PP and to the conclusion for a general 
discussion. 

4.2 Activated dynamic and scaling 

Another way to tackle and describe the aging dynamics that has been developed 
in the last twenty years is based on a growing correlation length, scaling and 
renormalization group ideas. The physical picture that one has in mind is similar 
to the one discussed for the aging dynamics of the Ising model. After a quench 
to the low temperature phase domains separating the two (or more) phases grow. 
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However, because of quench disorder, the domain walls have to overcome (free- 
energy barriers in order to grow. The scaling hypothesis is that at a timescale t 
the typical linear domain size is £ (t) and that correlation functions and responses 
at different times can be written in terms of a unique scaling function up to a 
length rescaling. For example, as for the non disordered Ising model the spin-spin 
correlation function is expected to read: 

<0(x,t)0(^O>=/(^^,|^) • (12) 

On this lengthscale, £(£), the typical (free-)energy barrier that domains have to 
overcome in order to grow is assumed to scale like As a consequence, using the 
Arrhenius law relating time and energy barrier one finds that £(t) oc (Tlogt) 1 ^, 
i.e. a very slow activated aging dynamics. Note that the previous expression is 
correct only in the regime 7^ 1. Also in this case one expects a pseudo- 
equilibrium regime, corresponding formally to t,t' >> 1 with |^ = 1 and 
physically to the equilibrium relaxation in the regions inside the domains. A 
systems in which this description is believed to be correct is the Random Field 
Ising model, RFIM, (ferromagnetic Ising model with independent random fields 
on each site). This system has a ferromagnetic phase transition for a not too high 
value of the variance of the random fields (their mean is zero). After a quench at 
low temperature domains separating the positive magnetization phase from the 
negative magnetization phase grows following the laws discussed above. This de- 
scription of the aging dynamics has been proposed for many different disordered 
systems like disordered ferromagnets, spin glasses, etc. ... In some cases, like 
for spin-glasses, this description of the aging dynamics is still matter of debate. 
Notice also that, at least superficially, it is very different from the mean-field one 
described in the previous section which is not based at all on activated dynamics. 
In the following I shall present a sketchy analysis of the off-equilibrium dynamics 
of the Id RFIM |7j. This model is an instructive example because the above 
heuristic predictions can be indeed obtained analytically and their physical con- 
tent can be easily grasped. The Hamiltonian is: 

2V-1 N 

H = - J ^2 S i S i+l - hiSi 
i=l i=l 

where hi are independent, say Gaussian, random variables with zero mean and 
variance g = hf. At zero temperature, the Imry-Ma argument implies that the 
region of + and — have a typical size Ljm = 4J 2 /g. The reason is simple: over 
a length / the energy gain obtained aligning the spins with the random fields is 
oc —y/gl whereas the energy loss for creating a domain wall between a region of + 
and — spins is 2 J. These two contributions match precisely at I = Ljm- Thus, we 
find that the RFIM in one dimension has of course no phase transition; however 
as long as the typical distance between domains during aging is less than Ljm 
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(or its finite-temperature counterpart) the system does not know it and the aging 
behavior is expected to be described by the activated dynamics scaling scenario 
described before. Note that this is also an example of interrupted aging in which 
if one waits long enough the system eventually equilibrates. 
The way to tackle the aging dynamics of the ID RFIM is to rewrite its Hamilto- 
nian in terms of positions of domain walls. Let me call a and b respectively the 
domains +|— and — 1 + , N a and N b their number, and and 6j the position of the 
ith domain of type a and b. Using this variables one can write the Hamiltonian 
as 

N N a N b 

H = -J(N -l)-J2h t + 2J(N a + N b ) + J2 V(oi) - £ V(h) 

i=l i i 

where V(x) = — 2^ =1 /jj. Furthermore, a Glauber dynamics on the spins (i.e. 
with a transition rate w(si — > — Sj) = e _/3AE / (e~^ AE + e +f3AE )) implies a reaction- 
diffusion dynamics for the domains in which a couple domains can be created on 
two neighboring sites with energy cost 4J ± hi, diffusion takes place with an 
energy cost ±2/ij and the domains annihilate on neighboring sites with an energy 
cost — 4 J ± hi. We are interested in the aging dynamics after a quench to a very 
low temperature T « J and with g « J 2 . This is the case in which L IM » I 
and the activated dynamics scaling picture should hold. In this case the effective 
dynamics for the domains is: a domains diffuse within the external potential V, 
b domains diffuse within the external potential —V and a and b annihilates when 
they meet. The activated dynamics scaling picture now becomes very concrete. 
At time t domains are at a typical distance £(£). This length is just the typical 
distance ds(t) on which a particle diffuses in the interval of time t in the random 
potential V (or —V, the sign does not matter). The main reason is that on 
lengthscales smaller than ds(t) there should be just one domain because if there 
were others then they would have meet and therefore they have been annihilated 
on timescale t. Instead on lengthscales larger than d s (t) all the domains that 
were present before the quench are still there, thus £(t) = ds(t). 
The problem of determining ds(t) and analyzing the motion before annihilation 
of a and b 'particles' is very well-known because is a discrete version of the Sinai 
Model: a particle diffusing with a Langevin dynamics in a Gaussian independent 
random force field. Indeed integrating the force to get the potential one finds 
precisely V(x) = Y%=i fi where fi are independent Gaussian random variables. 
It has been proved by Sinai that d,s{t) oc (Tlogt) 2 . There is also a simple 
argument to heuristically get that ip — 2: the potential V(x) in an interval of 
length d is like a Brownian motion (in which V is the position and x the time), 
thus the typical excursion of V around its boundary values is y/d. Hence, using 
the Arrhenius law we get precisely ds(t) oc (Tlogt) 2 . 

Finally, the other nice thing about the dynamics of the ID RFIM is that it can 
be really fully analyzed within a real space renormalization group procedure that 
provides an analytical derivation of the activated dynamics scaling picture and 
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of many other results This RG procedure has been introduced by Dasgupta 
and Ma and for disordered quantum spin chains and it has been developed and 
put on a firmer theoretical ground by D.S. Fisher and subsequently applied to 
many different systems [Hj. 

I will not continue further the presentation of this very appealing scenario for the 
aging dynamics and I refer to reviews for other details [HI E] and to the conclusion 
for a further discussion. 

4.3 Trap Model 

The trap model for aging dynamics has been introduced by Bouchaud in || and 
extensively developed since then pQ. The trap model is not a microscopic model. 
It is however very interesting because on the one hand it has been used and it 
is still used a lot on a phenomenological level to interpret results of experiments 
and simulation and, on the other hand, it provides a coarse grained description 
of different microscopic models. 

The simplest version jS] is defined as follows: there are N possible states plotted 
as wells in Fig. El The dynamical evolution is determined by stochastic jumps: 
each well acts as a trap for the system. The trapping time is determined by the 
Arrhenius law: r a = r exp(+(F — F a )/T) where F a is the free energy of the 
"state" (or well, trap) a and F is a threshold free-energy state, see Fig. El Once 

F 11 



F„ 



Figure 6: Schematic representation of the configuration space for the trap model, 
the system is escaped from a trap, it falls completely at random into another 
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trap. The rationale behind these rules is that generically for disordered systems 
there exists many different metastable states which, at least in some regimes, 
are surrounded by large energy barriers. The complete absence of geometry is of 
course the simpler assumption one can do. More complicated structures of the 
configuration space has been considered. The other essential ingredient of the 
model is the distribution of free-energies that is taken exponential with a certain 
parameter x: 



This choice is motivated by the presence of a similar distribution in mean-field 
disordered systems like the Random Energy Model (in which x = T/T 9 )pQ and, 
on general grounds, by the extreme value statistics for deep free-energy states [T|. 
The resulting distribution of trapping times ip(r) reads: 



The important feature of this law is that the average trapping time < r > is 
infinite for x < 1, i.e. T < T g if we take x as proportional to the temperature. In 
this regime the system, taken at time zero uniformally distributed over the traps, 
never reaches the equilibrium and ages forever. In this case the nature of the 
aging dynamics is particularly clear 0: the time t w after n jumps equals Yh=i T i 
where 7$ are independent variables with probability law ip(r). It is a well-known 
result that since the first moment is infinite the central limit theorem does not 
apply, t w is distributed with a Levy Law and its typical value is of the order of the 
maximum trapping time encountered during the evolution. At a given waiting 
time t w the sum n + T2 + . . . is dominated by its largest term which is of the 
same order of t w . In this case it is very clear that the system ages because the 
relaxation time is set by the age of the system itself. 

I will not present further results on the trap model and I refer to [§1 E] for a 
wider presentation. I would like just to stress that the Trap Model provides a 
coarse grained description of the Sinai model in a biased potential jS] as well as 
the aging dynamics (on extremely long times, diverging with the system sizes) 
of the Random Energy Model fTO]. These results suggest that it is not only a 
phenomenological model and it might play in the future an important role in the 
description of the finite dimensional extension of the mean-field aging. 

5 Outstanding issues 

In this section I would like to point out what are in my mind some outstanding 
theoretical issues in aging of glassy systems. It will be certainly amazing to look 
at this list in a couple of years to see how things turned out to be very different. 
Beyond the mean-field theory of Aging. 

Length-scales, Landau Theory and fluctuations. As for critical phenomena one 



P(F) = x/Texp(x{F - F )/T). 
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can identify three steps that have to be followed to analyze the physics of finite 
dimensional systems. First, one can solve an infinite dimensional system. This 
has been done for aging and these studies indeed lead to a mean-field theory of 
aging 1J. A second important step is set up a Landau theory. Although the 
results still are mean-field like (i.e. no fluctuations are taken into account) this 
provides hints of why results might be universal (at least in higher enough di- 
mension) of which type of correlation lengths are diverging (as we discussed in 
the introduction if the aging is not interrupted one expect a diverging length) 
and it helps to clarify what are the important fluctuations that one has to take 
into account. I think that it is a very important point that will be studied in 
detail in the future, some results have been already obtained in [TT] . 
Beyond mean-field theory: activated processes. For mean-field systems conjec- 
tured to be related to structural glasses the thermodynamic as well as the equi- 
librium and off-equilibrium dynamics are expected to change drastically in finite 
dimension. The main reasons is that the physics of these models is dominated 
by the existence of an infinite number of thermodynamic states and, instead, it 
is clear that in finite dimension this exponential proliferation of thermodynamic 
states cannot take place. How this is going to change the physics (statics and, 
specially, dynamics)? How this translate in terms of real space structure? What 
remains of the off-equilibrium dynamics described in terms of energy landscape 
picture? These are certainly very important and pressing questions that, hope- 
fully, will receive a lot of attention in the near future. 
Renormalization group for the activated dynamic scaling. 
Although in one dimension there are examples of activated dynamic scaling that 
can be derived and fully analyzed analytically within a RG procedure the ex- 
tension to higher dimensions remains an important open problem and certainly 
a major challenge for the future. Recent interesting results have been already 
obtained in [121 ITS] . 

Generalization of the Fluctuation-Dissipation relation in finite dimen- 
sion 

This question has already attracted a lot of attention. In particular one would 
like to know if and how the FDR is generalized in finite dimension in comparison 
with mean- field predictions and if the resulting T e ff can be really interpreted as 
an effective temperature. Furthermore also the generalization of FDR to driven 
off-equilibrium systems as for example granular media is an interesting related 
question. Despite many simulations and some experiments have been done these 
questions are not completely settled and certainly they will continue to give rise 
to many works in the future. 

Aging in glassy systems without quenched disorder 

There are systems like structural glasses that have a very interesting slow dy- 
namics that is still a puzzle despite many years of theoretical and experimental 
investigations. The understanding of their off-equilibrium dynamics is therefore 
a challenge for the future. 
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From this point of view it is very likely that in the future much more attention will 
be devoted to the off-equilibrium dynamics of simple statistical mechanics model 
of glassy systems as Kinetically Constrained Models, Lattice Glasses, . . . (see [H] 
for a review). 

Quantum off-equilibrium dynamics What about quantum glassiness? This 
is certainly a very interesting subject that will probably be the target of many 
investigations in the future. To me very interesting questions are: what are the 
properties of off-equilibrium dynamic (aging) close to a quantum critical point? 
What are the properties of the aging behavior induced by a quench (that goes 
through a quantum phase transition) due to a sudden change in the control 
parameter governing quantum fluctuations? What is a quantum glass without 
disorder? 

Some of these questions have been recently addressed in [Tol ITS]. 

I will stop here this list that could go on much further. I think I have described 

already enough challenging problems. 

6 Suggested readings 

Suggested readings that are useful to have a deeper understanding of the aging 
dynamics and the theoretical ideas developed to understand it are: 
J. Kurchan and L. Laloux, J. Phys. A 29, 1929 (1996) for the energy descent 
picture of aging. 

L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993) for the aging 
dynamics in mean field glassy systems. 

D. S. Fisher and D. A. Huse, Phys. Rev. B 38, 373 (1988), Phys. Rev. B 35, 
6841 (1987) for the activated dynamic scaling. 

J. -P. Bouchaud, J. Phys. I France 2, 1705 (1992) for the Trap model. 
This list is clearly not exhaustive. 
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